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ABSTRACT 

Unusual patterns of recombination emission from gas ionized by metal-free stars may distinguish early star- 
forming galaxies from their present-day counterparts. This pattern arises from the harder ionizing spectrum ex- 
pected from metal-free stars, which strongly enhances the strength of He II recombination lines. Our calculations 
indicate that line fluxes of He II A 1640 and A4686 are sufficiently large to be detected by narrowband and spectro- 
scopic searches for high-redshift emission-line sources at z ~ 5 using current instruments. An unknown fraction 
of Lya emitters may harbor low-metallicity or metal-free stars. As the predicted He II A1640 flux is comparable 
to and may exceed hydrogen Lya, searches for high-redshift galaxies should consider He II recombination lines 
as possible identifications for single emission lines in observed spectra. Spectra of metal-free stars may show 
both H I and He II emission lines, improving the constraints on their redshift and identification. We assess the 
considerable uncertainties that affect our expectations for the detection and identification of true first-generation 
stars with present search techniques, including the role of stellar mass loss in spectral evolution and the confusion 
of ionization by primordial stellar sources and AGN in the early universe. 

Subject headings: stars: early type — intergalactic medium — cosmology: theory 



1. INTRODUCTION 

The first stars have long attracted interest from a cosmolog- 
ical point of view, and recent detailed studies of their prop- 
erties have renewed interest in their cosmological importance. 
Within the constraints of big-bang nucleosynthesis, these stars 
must have been metal free, which modifies their structure in 
fundamental ways: metal free stars have smaller radii, hotter 
cores, and higher effective temperatures than their counterparts 
of equal mass but finite metallicity (Tumlinson & Shull 2000). 
These stars are leading candidates for the sources that reionized 
the intergalactic medium (Gnedin 2000; Haiman & Loeb 1997), 
and they produced the first heavy elements that enriched sub- 
sequent stellar generations. Tumlinson & Shull (2000) found 
that these stars, here called Population III or simply metal-free, 
have harder spectra and emit more ionizing photons in the He II 
ionizing continuum than do stars of typical solar or sub-solar 
metallicity. Thus, they may dominate He II reionization in the 
high-redshift universe. 

Distinctive spectral features of their surrounding nebulae 
may provide a means of detecting the first stars, during and af- 
ter the epoch of reionization. Tumlinson & Shull (2000) noted 
that the hard spectra of Pop III stars produce large He III re- 
gions, which may emit detectable He II recombination emis- 
sion. Recombination lines of He II at A 1640 (n = 3 — ► 2), A3203 
(n = 5 — > 3), and A4686 (n = 4 — > 3) are particularly attractive 
for this purpose, because they suffer minimal effects of scat- 
tering by gas and decreasing attenuation by intervening dust. 
Thus, an assessment of the near-term prospects for discover- 
ing metal-free stellar sources is timely. However, uncertainties 
about the primeval initial mass function (IMF), the importance 
of stellar mass loss at zero metallicity, and the features of stellar 
evolution must be addressed. 

In this Letter, we consider the possibility of detecting metal- 
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free stars using He II recombination emission, and the uncer- 
tainties inherent in this technique. Main-sequence mass loss, 
not well understood from an evolutionary standpoint, can af- 
fect the detectability of the first stars via nebular emission lines. 
Emission line ratios and line profiles may discriminate between 
stars and active galactic nuclei (AGN) as sources of ionization. 
We address these issues in the context of existing tracks and 
assess their importance to the detectability of metal-free stars. 

In § 2 we derive an expression for the observed He II 
emission-line flux of a metal-free stellar population in terms 
of the star formation rate and evolutionary properties. In § 3 we 
use evolutionary tracks and model atmospheres to assess the ef- 
fects of stellar evolution on the detectability of metal-free stel- 
lar populations. In § 4 we compare our predictions to current 
observational searches for high-redshift stellar populations, in- 
cluding present-day emission-line techniques, and we discuss 
the ambiguity between emission lines produced by metal-free 
stars and AGN. In § 5 we summarize our results and comment 
on directions for future work. 

2. EMISSION LINES FROM METAL-FREE STARS 

We begin by estimating a relationship for the observed He II 
emission-line flux from metal-free stars as a function of star 
formation rate (SFR) and stellar evolutionary parameters. Us- 
ing the zero-age structure and atmosphere models presented by 
Tumlinson & Shull (2000), we use recombination theory to es- 
timate the total luminosity in the He II emission lines. We as- 
sume that few He II ionizing photons escape the galaxies where 
they are produced (f esc ~ 0) and that the sources themselves 
contain no dust. We use the Kennicutt (1983) law to relate the 
luminosity of the Ha line to the star formation rate, L(Ha) = 
(1.12 x 10 41 erg s _1 ) (SFR/M Q yr _1 ). For case B recombina- 
tion at T = 20,000 K (corresponding to a higher nebular tem- 
perature in low-metallicity gas) and «Hein/"Hii = 0.0789, we 
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find jAm/ha = 0.33, jmo/ha = 2.3, and him/ ha = 0.14 
(Seaton 1978). Kennicutt (1983) assumed a Salpeter-like IMF 
and used stellar evolution tracks for Population I massive stars. 
We assume a similar IMF and include a factor, / evol (of order 
unity), designed to account for the time evolution of the stel- 
lar ionizing continuum radiation, and defined to be unity if the 
evolution in He II ionizing photon rate for metal-free stars is 
identical to the evolution of H I ionizing photons used by Ken- 
nicutt (1983). The factor / evo ] can be evaluated using Pop III 
evolutionary tracks coupled with model atmospheres (see § 3). 
Using the results of Tumlinson & Shull (2000), we scale the 
He II ionizing photon production of zero-metallicity stars to the 
H I ionizing photon production implicit in the Kennicutt (1983) 
relation. We find that Pop III stars produce 10" 1 1 as many He II 
ionizing photons as the Pop I stars produce H I ionizing pho- 
tons. Scaling the luminosity of He II A1640 to Ha, we find: 



L 1640 = (4.2 xl0 40 ergs- 1 )/ evol 



SFR 

M yr- 



(1) 



For deceleration parameter qo = 0.5, the flux of this source is: 

(4.1 x 10" 18 erg cm" 2 s" 1 )^, f evol SFR 
F mo (z) = B - , (2) 



[(l+z)-(l+z)5] 2 

where is the Hubble constant in units 65 km s" 1 Mpc" 1 . The 
flux of He II A4686 is 7.1 times lower than Fi64o, and the flux 
of A3203 is 16.2 times lower. (Hereafter we concentrate on the 
two stronger lines.) This relation is plotted in Figure 1 for SFR 
= 5 and 20 Mq yr" 1 and / evo i = 0.4 and 2.0; these choices of 
/evoi are justified in § 3. We also plot detection limits of three 
recent emission-line searches for high-redshift galaxies. Figure 
1 shows that, even for conservative estimates of / eV oi> metal-free 
stellar populations may be detectable to z = 2 - 5 for a reason- 
able range of SFR. If a detection can be made and its redshift 
accurately measured, the star formation rate can be constrained 
if /evoi is carefully calibrated by stellar evolution models. 

3. STELLAR EVOLUTION EFFECTS 

Stellar evolution at zero metallicity is uncertain and may en- 
hance or diminish the detectability of metal-free stars. Tumlin- 
son & Shull (2000) estimated the ionizing photon production 
from zero-age main sequence (ZAMS) metal-free stars of mass 
2-90 Mq using static stellar structure models and NLTE model 
atmospheres. The He II ionization produced by these stars is 
a direct result of their high effective temperatures. However, 
published evolutionary tracks of metal-free stars (Castellani et 
al. 1983; Chieffi & Tornambe 1984) show that these stars may 
evolve to cooler temperatures and larger radii over their life- 
times if they do not experience mass loss. Because mass loss 
may play a significant role in the spectral evolution of the star, 
it produces uncertainty in the interpretation of emission-line di- 
agnostics. 

With model atmospheres similar to those presented by Tum- 
linson & Shull (2000), Hubeny, Lanz, & Heap (2000) argue 
that line-blanketed, radiation pressure-driven winds are not ini- 
tiated for stars with Z <; 0.001, owing to the relative lack of 
metal line-blanketing in their atmospheres. Kudritzki (2000) 
draws similar conclusions. However, these groups have not, as 
yet, extended their models to the range T e ff > 60000 K occupied 
by metal-free massive stars. 

Because these stars radiate near the Eddington limit, they 
may drive winds with electron-scattering opacity. El Eid et al. 



(1983) used an empirical mass-loss prescription from Chiosi 
(1981) that scaled M to the Eddington luminosity with no ex- 
plicit dependence on metallicity. They found typical mass loss 
rates of M ~ 10" 5 M© yr" 1 for stars of 80 - 500 M , with lit- 
tle dependence on luminosity. At these rates, a star with initial 
mass 80 Mq will fall to 40 Mq by the end of hydrogen burning, 
when it begins an excursion to hotter T e ff in the HR diagram. 
As its outer layers are lost, hotter regions below are exposed, 
and the interior structure adjusts to the changing mass and pho- 
tospheric conditions. At late times, T e s increases by 0.2 dex 
over their ZAMS, increasing the production of He II ionizing 
photons by an order of magnitude. 

To evaluate the importance of stellar evolution and mass loss 
on the production of He II ionizing photons, we make two lim- 
iting assumptions: (1) mass loss from Population III stars is 
negligible, with no effect on their evolution; (2) mass loss is 
important and affects the structural and spectral evolution of 
the star. We use the 25 Mq track calculated by Castellani et 
al. (1983) to represent the first case, and the 80 Mq mass-loss 
track by El Eid et al. (1983) for the second. These tracks are 
not intended to replace the stellar models on which our flux es- 
timates are based. Rather, they represent extremes of mass loss 
used to constrain the parameter space of / evo i- If mass loss is 
unimportant to the evolution of metal-free stars, then the gen- 
eral trend to lower T eB seen in the 25 Mq track should hold for 
stars at all mass and will favor low values of / evo i- Conversely, 
if mass loss substantially affects the later evolutionary phases of 
metal-free stars, then the trends toward higher T e ff at late time 
will favor larger / evo i- 

Model atmospheres were calculated with the TLUSTY code 
(Hubeny & Lanz 1995) for a set of points along the evolutionary 
tracks, placed to capture the time evolution of T eB and surface 
gravity g. We define a time-averaged flux parameter, 



/e,= 



G,(0)f PSN ' 



(3) 



where 2, = o,L2 are the ionizing photon production rates for H I, 
He I, and He II, respectively, and f P sN is the pre-supernova life- 
time of the star. The 25 Mq star evolves to lower T e a during 
H burning and experiences a burst of He II ionizing photons 
near the onset of He burning. The 80 Mq star loses over half 
of its mass during its lifetime and makes two brief excursions 
to high T e ff and high Q2. For the constant-mass (25 Mq) track, 
we derive / Qo = 1.3, / Ql = 1.1, and f Ql = 0.32. For the 80 M 
track with mass loss, we find /q = 1.2, /q, = 1.1, and /q 2 = 
1.80. For comparison, we estimate that /q = 0.7, /q, = 0.4, and 
/q 2 = 0.3 are characteristic of the Chiosi, Nasi, & Sreenivasan 
(1978) Population I tracks used by Kennicutt (1983). Because 
/evoi compares the evolutionary trends in the He II ionizing flux 
from metal-free stars with the time average in H I ionizing flux 
from Pop I, our estimate is given by / evo i = /q 2 (Pop III)//q 
(Pop I). Thus, on timescales long compared to the evolution of 
massive stars, / evo i may range from 0.4 to 2.6. The two limiting 
cases imply different behavior on short timescales. For M = 0, 
the He II emission fades and H I emission brightens over time 
as the stars evolve to cooler temperatures. For M =/0 the He II 
emission from a Pop III cluster increases over time, as stars 
of successively lower mass exhaust their central hydrogen and 
make excursions to hotter T e g. 

Our calculation illustrates a key feature of the ionizing pho- 
ton production by the first stars: emission-line diagnostics 
evolve with their stellar populations, and this must be con- 
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FIG. 1 . — Flux of He II A1640 line as a function of source redshift (eq. 2). The dashed and solid curves correspond to / evo i = 0.4 and 2.0, respectively. The 
upper and lower curves in each pair correspond to star formation rates 20 and 5 Mq yr 1 . The flux for A4686 is 7.1 times lower than F 1640 . The horizontal solid, 
dotted, dash-dot, and dashed lines show the limits of current Lya emission-line surveys, as described in the text. The flux limits have been corrected for the energy 
difference between 1216 A and 1640 A. At the top, we mark the ranges of redshift probed by the He II A 1640 (above) and A4686 (below) lines for optical and 
infrared searches. Optical searches for z = 1 -5 are feasible now, and infrared searches may be possible in the near future. 



sidered in the planning and interpretation of observational 
searches. A single snapshot of a metal-free stellar population 
in its He II recombination lines would not distinguish between 
the M = and M ^ cases, because of the degeneracy in the 
parameters / evo i and SFR (eq. 2). Measuring the star formation 
rate of metal-free stars with He II emission requires that /q 2 be 
carefully calibrated by stellar evolution calculations. 

The evolution of the ionizing photon production may pro- 
vide a means of distinguishing between the two cases presented 
above. The time evolution of g; for the limiting cases of M = 
and M / illustrates a key difference between the two possi- 
bilities. The 25 Mq star at constant mass exhibits most effi- 
cient He II ionization at the beginning of its H-burning main- 
sequence, giving Q2/Q0 = 10" 4 . The 80 Mq star with M=/0, by 
contrast, achieves the maximum value Q2/Qo = 0.40 at the end 
of its life. This ratio is larger than the range of Q2/QQ = 0.05 - 
0.12 for the most massive stars studied by Tumlinson & Shull 
(2000) and Bromm, Kudritzki, & Loeb (2000). Thus, the two 
extreme cases can be distinguished by observed He II to H I 
ionization ratios that exceed those achieved by constant-mass 
stars. 

4. COMPARISON WITH OTHER EMISSION-LINE DIAGNOSTICS 

The relatively bright He II emission-line fluxes predicted for 
metal-free stars raise the possibility of detection with present- 
day broadband and spectroscopic searches. In particular, we 
examine the possibility that He II A 1640 emitters have already 
been observed and mistaken for Lya emitters. Recent searches 
for high-redshift galaxies use narrow-band imaging and spec- 
troscopic techniques. These emission lines are believed to be 
Lya because deep spectra reveal no other galactic emission 
lines ([O III] A5007, [O II] A3727, Ha), and these identifica- 
tions for the single line can be eliminated by the absence of 



corresponding emission. 

In Figure 1 we plot the detection limits of three current Lya 
searches for high-redshift galaxies. Hu, Cowie, and McMa- 
hon (1999) quote a 5cr detection limit 1.5 x 10" 17 erg cm" 2 s" 1 
for narrowband imaging and 5.0 x 10~ 18 erg cm" 2 s" 1 (5a) for 
spectroscopic searches. In a 1.5 hour Keck exposure, Stern et 
al. (2000) achieved a limiting flux 6.0 x 10" 18 erg cm" 2 s" 1 in 
a 1" slit for an unresolved line. Rhoads et al. (2000) achieved 
detections ranging between 1.8-2.6 x 10" 17 erg cm" 2 s" 1 us- 
ing combinations of imaging and spectroscopy. As shown in 
eq. 2, such a detection of He II lines would measure the prod- 
uct of SFR and / evo i- Therefore, careful calibration of / evo i with 
stellar evolution models could break the degeneracy and inde- 
pendently constrain the star formation rate. 

Hu, McMahon, & Cowie (1999) reported the detection of a 
single line in the spectrum of a luminous galaxy, using a combi- 
nation of narrowband imaging and spectroscopic observations. 
Identifying this line as Lya, they infer z = 5.74, a total line flux 
1 .7 x 10~ 17 erg cm" 2 s" 1 at 8 195 A, and a star formation rate of 
19 M Q yr" 1 . If this line is He II A1640 at z = 4.0, the implied 
SFR would be 80, 32, and 16 Mq yr" 1 for / evo i(G2) = 0.4, 1.0, 
and 2.0, respectively. Typical Lyman-break galaxies have in- 
ferred SFR ~ 10- 100 Mq yr" 1 (Steidel et al. 1999). Thus, the 
SFR by itself does not rule out identification as He II A 1640 for 
implied values of /evoi consistent with the M^0 case. The im- 
plied total flux in A4686 is 2.3 x 10" 18 erg cm" 2 s" 1 at 2.34 fim, 
below the K-band sensitivity limits of the Hu et al. (1999) sur- 
vey. The flux predictions in § 2 suggest that this line might be 
He II A1640 at z = 4.0, were it not for the presence of a spectral 
break across the line produced by intergalactic Lya absorption. 
The absence of the break may provide key spectroscopic evi- 
dence for identifications of He II recombination lines. Without 
a break, the U, B, and V dropout techniques commonly used 
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to discover high-redshift galaxies (Steidel et al. 1999) may not 
select these objects for spectroscopic followup. Because Qq for 
Pop III is similar to that of Pop I (Tumlinson & Shull 2000), the 
observed fluxes of He II A 1640 and H I Lya should be compa- 
rable. The predicted line fluxes and the separation of Lya and 
He II A 1640 by 424(1 +z) A suggests that some detections may 
show two emission lines, confirming identification and redshift 
in the range z = 2-5. Future spaceborne missions such as NGST 
may also be used to detect high-redshift AGN and metal-free 
stars via associated He II A 1640 and A4686 nebular emission 
(Oh, Haiman, & Rees 2001). 

It remains an issue, however, how the detection of a single 
nebular line would permit one to discriminate between stellar 
and AGN sources of ionizing radiation. Intrinsic line widths 
may not be sufficient to resolve the confusion. Recent surveys 
of Wolf-Rayet galaxies (Schaerer et al. 1999 and references 
therein) show broad He II A4686 emission from stellar winds 
and narrow nebular emission. These lines are usually detected 
at fluxes of 1 - 4% of the observed H/3 fluxes from these galax- 
ies, and the broad wind features are generally brighter than the 
nebular emission. Spectroscopic Lya surveys (Hu et al. 1999) 
typically use low resolution, R ~ 1000, to achieve their required 
sensitivity. Identifying stellar and AGN ionization may require 
resolution of ~ 100 km s -1 to distinguish broad, cuspy line pro- 
files of AGN lines from narrow nebular emission due to massive 
stars. The He II to H I line ratios for metal-free stars may pro- 
vide a further discriminant. The M ^0 case (§ 3) gives a maxi- 
mum value of Q2/Qo = 0.40. For a power-law AGN spectrum, 
/ ~ v" Q , this ratio is Q2/Q0 = 4" Q . For the composite radio- 
quiet AGN spectrum (a = 1.8) from Zheng et al. (1997), this 
ratio is Q2/Q0 = 0.08. However, distinguishing metal-free stars 
from AGN with harder spectra than the power-law composite 
may be more difficult. If mass loss is important to the evo- 
lution of metal-free stars, this maximum ratio technique may 
provide a convenient means of distinguishing them from AGN 
in the early universe. The prospects for distinguishing stellar 
and AGN ionization improve for cases where lines of both H I 
and He II are observed. 

As a potential discriminant between AGN and metal-free 
stars, IGM ionization ratios of He II to H I may show dif- 
ferent signatures of AGN or stellar ionization. However, the 
hard AGN-like ionizing continuum of metal-free stars may con- 
fuse the issue. Fardal, Giroux, & Shull (1998) define the col- 
umn density ratio 77 = N^ii/Nhi- Ultraviolet studies of the 
He II Gunn-Peterson absorption toward quasars at z = 2.7-3.3 
(Davidsen, Kriss, & Zheng 1996; Reimers et al. 1997; Hogan, 



Anderson, & Rugers 1997) suggest that the IGM reionizes in 
He II at z ~ 3. The ionizing spectrum required to match the 
optical depths, thi and THeii, at z < 3 is relatively soft, with 
r\ = 100-200, in contrast to the harder spectra from unfiltered 
AGN and metal-free stars, which yield 77 w 30. If metal-free 
stars are implicated in the He II reionization, then observational 
searches for the characteristic He II emission lines (A1640 and 
A4686) from starburst galaxies at z = 3 - 5 would be extremely 
helpful. These lines should be measured together with high- 
ionization metal absorption lines (C IV, S IV), to assess any 
shift in the quality of the ionizing spectrum at z ~ 3 (Songaila 
& Cowie 1996). Discriminating stellar from AGN sources may 
depend on detailed line ratio predictions from accretion models 
(Hubeny et al. 2000) and stellar population synthesis. The max- 
imum line ratio technique sketched above must be extended to 
account for IMF and evolution effects, but it holds promise if 
complete H I and He II recombination spectra are observed. 

5. SUMMARY 

In summary, our main conclusions are the following: 

1. Predicted fluxes for He II emission lines from nebulae 
ionized by metal-free stars are comparable to the de- 
tection limits of present-day narrowband and spectro- 
scopic searches for high-redshift galaxies. Observational 
searches for these stars at z ~ 3 - 5 are feasible now at 
optical wavelengths. Higher redshift surveys may need 
to await space-based infrared instruments. 

2. Mass loss by metal-free stars may restrict or enhance 
their detectability with emission-line techniques. Careful 
calibration with detailed evolutionary tracks and theoreti- 
cal work on mass loss from metal-free stars are necessary 
to refine this program. 

3. If multiple line detections are made for a single object, 
He II and H I emission lines may provide important di- 
agnostics of the stellar parameters. These line ratios 
may provide a means of distinguishing stellar and AGN 
sources even if the line profiles are unresolved. 
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